
Europaisches 
Patentamt 



10 / 5 0 4 9 l3..-'MEP .P 3 / 0 1 9 4 5 



European 
Patent Office 



Office europeen 
des brevets 




Bescheinigung Certificate 



Attestation 



Die angehefteten Unterla- 
gen stimmen mit der 
ursprQnglich elngereichten 
Fassung der auf dem nSch- 
sten Biatt bezeichneten 
europ§ischen Patentanmel- 
dung aberein. 



The attached documents Les documents fix^s d 
are exact copies of the cette attestation sont 
European patent application confbrmes d la version 
described on the following initialement d6posee de 
page, as originaily filed. la demande de brevet 

europeen sp6clfi6e d la 
page suivante. 



Patentanmeldung Nr. Patent application No. Demande de brevet n"* 

02388019.8 



I PRIORITY 
, DOCUMENT 

I SUBMITTED OR TRANSMITTED IN 
' COMPLIANCE WITH RULE 17.1(a) OR (b) 



Der Prdsident des Europdischen Patentamts; 
Im Auftrag 

For the President of the European Patent OfRce 

Le Pr&ident de I'Office europeen des brevets 
p.o. 



R c van Dijk 



BEST A\/A».ARI F COP*^^ 




# 



Europaisches 
Patentamt 



Euroj 
Patei 




f Office europten 

Ice des brevets 



PCT/EP 0 3/01945 



Anmeldung Nr: 
Application no.: 
Demande no: 



02388019.8 



Anme1 detag: 
Date of filing: 
Date de d^pfit: 



07.03.02 



Anmel der/Appl 1 cant( s )/Demandeur( s ) : 
Telefonaktiebolaget L M Ericsson (Publ) 



126 25 Stockholm 
SUEDE 



Bezelchnung der Erf 1 ndung/Tltle of the Inventlon/Ti tre de 1' invention: 
(Falls die Bezelchnung der Erflndung nicht angegeben 1st. slehe Beschreibung. 
If no title Is shown please refer to the description, 
SI aucun titre n'est 1nd1qu€ se referer & la description.) 

Soft value calculation for multilevel signals 

In Anspruch genommene Prloriat(en) / Prlorlty(les) claimed /Pr1or1t6(s) 

revend1qu£e(s) ^ ^ 

Staat/Tag/Aktenzelchen/State/Date/Flle no./Pays/Date/Numfiro de d6p6t: 



Internationale Paten tkl as si f1 Rati on/Inter national Patent Classification/ 
Classification Internationale des brevets: 

H04L27/00 

Am Anmel detag benannte Vertragstaaten/Contracting states designated at da 
flllng/Etats contractants designees lors du dSpdt: 



AT BE CH CY DE DK ES FI FR GB GR IE IT LI LU MC NL FT SE TR 



02388019.8 

EPVEPO/OEB Form 1014. E - 01.2000 7001014 



2 




- 1 - 



Soft value calculation for multilevel signals 

This invention relates to digital communications systems and, more 
particularly, the generation of soft reliability values for multilevel signals. 

6 

Within the field of digital communications, multilevel modulation Is used to 
map a number of bit sequences to a signal alphabet comprising a number of 
signal symbols, i.e. a number of points in signal space. For example, a bit 
sequence may be mapped onto a point In a complex signal space. A signal 

10 alphabet of size M allows log2(M) bits to be mapped to each symbol. 
However, when symbols are received at a receiver, they may be affected by 
noise, thereby affecting the decoding of the signal when retrieving the 
transmitted bit sequence. If multilevel modulation is used in conjunction with 
channel coding, many channel decoders, such as iterative decoders based 

1 5 on the BC JR algorithm, require so-called soft bit values as an input A soft bit 
value con-esponds to a reliability value of a single bit being 0 or 1. 

Examples of multilevel modulation include multi-amplitude level modulation In 
Pulse Amplitude Modulation (JPm), multl phase level modulation in Phase 
20 Shift Keying (PSK), multl signal point modulation in Quadrature Amplitude 
Modulation (QAM). 

For example, an emerging technology for wideband digital radio 
communications of Internet, multimedia, video and other capacity-demanding 

25 applications In connection with the third generation of mobile telephone 
systems is the evolving Wideband Code Division Multiple Access (WCDMA) 
specified as part of the 3GPP standardisation organisation. Within this 
technology, High Speed Downlink Packet Access (HSDPA) Is provided 
Including a high speed downlink shared channel (HS-DSCH) which uses 

30 16QAM. In 16QAM for example, M=16. I.e. each symbol in the signal 
alphabet represents 4 bits. Future releases may comprise even larger 
constellation sizes such as 64 QAM. 

It Is known how to convert signal symbols to soft bit values by calculating all 
35 distances in signal space between the received symbol and all signal points 
of the signal alphabet. In particular, in order to obtain optimal perfomiance, a 
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likelihood ratio is calculated depending on corresponding sums of 
probabilities where the probabilities are functions of the calculated distances. 
It is further known that in the calculation of a likelihood ratio the sums of 
probabilities may be approximated by the dominant contributions to the sums 
5 of probabilities In the likelihood ratio (A. J. Vlterbl, "An intuitive justification and 
a simplified Implementation of the MAP decoder for conventional codes". 
IEEE Journal on selected areas In communications, 16(2), February 1998). 

Even though this approximation significantly reduces the computational 
1 0 complexity while only causing a negligible loss In performance, It is a problem 
of the above prior art method that It still requires a calculation of all distances 
to all the signal points In order to determine which two are actually needed for 
the calculation of the likelihood ratio. For example, In a 16QAM modulation. 
16 distances have to be calculated for each received symbol. In particular, if 
15 a high rate of symbols need to be decoded, e.g. several hundred symbols per 
millisecond, the above problem Is particularly severe. 

The above and other problems are solved when a method of generating a 
reliability value for a received multilevel signal In relation to a number of 
20 predetermined signal symbols each associated witii a corresponding bit 
sequence including a first bit position; the reliability value being Indicative of 
likelihood information of receiving said multilevel signal is characterized in 
that the method comprises the steps of 

- identifying a first one of the number of signal symbols as being closest to 
25 the received multilevel signal; and 

- estimating the reliability value based on stored Information related to a 
second one of the number of signal symbols being closest to the first 
signal symbol and corresponding to a different binary value at the first bit 
position of the respective associated bit sequence than the first signal 

30 symbol. 

Consequently, according to the invention, the closest signal point needed for 
calculating the likelihood ratio is determined first, such that only the 
corresponding distances to the Identified signal points need to be 
35 determined. In this way, the computational complexify is reduced 
significantly, since not more than two distances need to be calculated In 
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order to determine the likelihood ratio for each bit. Hence, the at>ove method 
is particularly well-suited for low-complexity implementations in mobile 
receivers, as It reduces the required computational resources. 

According to a prefen-ed embodiment of the invention, the stored information 
comprises an identification of the second signal symbol, and the step of 
estimating the reliability value further comprises the steps of 

- determining a first distance between the received signal and the first 
signal symbol; and 

- determining a second distance between the received signal and the 
second signal symbol. 

Hence, given a signal constellation, as for each signal symbol and for each 
bit position it is known which one of the other signal symbols having an 
16 opposite value at that bit position is closest to that signal symbol, this 
information may be stored. Therefore, by a simple look-up operation the 
closest signal symbol may be identified and the distance between the 
received signal and the closest signal symbol may be calculated. 

20 According to another prefen-ed embodiment of the invention, the stored 
infomiation comprises a second distance between the first signal symbol and 
the second signal symbol, and the step of estimating the reliability value 
further comprises the step of detenmlning a first distance between the 
received signal and the first signal symbol. Hence, alternatively or 

25 additionally to storing an identification of the signal symbol with opposite bit 
value at a given bit position being closest to the first signal symbol, the actual 
distance between the first and the second symbol may be pre-calculated and 
stored. Once the first signal symbol is identified, this distance may be looked 
up and used as an approximation for the distance between the received 

30 signal and the second signal symbol. Consequently, a further reduction in 
computational complexity is achieved, since only one distance has to be 
calculated. 

When the step of estimating the reliability value comprises the step of 
35 detemiining a polynomial function of the first distance and the second 
distance, multiplied by a predetermined constant, the computational 
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complexity Is further reduced, as no logarithm needs to be calculated. In one 
embodiment, the polynomial function is a difference of the squared distances. 

According to yet another prefenred embodiment of the invention, the stored 
5 information is indicative of a number of functional relations between the 
received multilevel signal and the reliability value, and the step of estimating 
the reliability value further comprises the step of selecting a functional 
relation of said number of functional relations dependant on the first signal 
symbol and the first bit position. Hence, the calculation of the likelihood value 
10 only comprises the step of calculating the corresponding stored function. 
Preferably, the functional relationship Is a linear function of a signal 
component, thereby reducing the calculation to a multiplication operation and 
an adding operation. 

15 According to yet another preferred embodiment, the stored Infonnatlon 
comprises, for each signal symbol and bit position an approximation of the 
corresponding reliability value. Hence, according to this embodiment, an 
approximation of the reliability value may be directly looked up once the 
closest signal symbol Is identified, thereby providing a computationally very 

20 efficient method which eliminates the need of online distance calculations. 

Preferably, the stored infomnation is stored In a look-up table indexed by the 
number of signal symbols and the bit positions of the number of bit 
sequences, thereby providing fast access to the Information. 

26 

In one embodiment of the invention, the method further comprises the step of 
providing the rellabiiify value as an input to a decoder, e.g. an iterative 
decoder using the BCJR algorithm or any other decoder using soft values as 
an input. It is an advantage of the invention that it provides an accurate and 
30 resource-efficient approximation of soft values as an input to such decoders. 

The first signal symbol may be identified by comparing the signal 
components with predetermined thresholds or decision boundaries, for 
instance by means of a slicer, i.e. a circuit which compares a signal with 
35 predetermined thresholds. Hence, a fast and computationally inexpensive 
method is provided for identifying the closest signal symbol without the 
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necessity of calculating all distances between the received value and all 
signal symbols. 

It is a further advantage of the invention that cost-effective, standard 
5 components may be employed when implementing a method according to 
the Invention. 

When the likelihood information comprises a log-likelihood ratio, a high 
performance quality is achieved, as the use of a log-likelihood corresponds a 
10 theoretically optimal way of calculating soft reliability values. However, other 
methods of calculating likelihood Infomiation may be employed, such as a 
log-llkellhood of the signal power. 

In a preferred embodiment of the invention the step of identifying the first 
1 5 signal symbol as being closest to the received multilevel signal comprises the 
step of identifying the first signal symbol as being closest to the received 
multilevel signal with respect to a Euclidean distance measure In a signal 
space, as the Euclidean distance Is directly related to the probabilities of a 
likelihood calculation. Alternatively, other suitable known metrics may be 
20 used instead Euclidean distances. 

The signal space may be a real or complex signal space. For example, in 
QAM modulation two amplitude-modulated signals are transmitted on a 
single can-ier. but shifted In phase by 90 degrees. Hence, the resulting signal 
25 points may be represented in the complex plane representing the so-called 
in-phase (I) and quadrature (Q) components of the QAM signal. 

When the number of signal symbols are associated with the number of bit 
sequences such that the bit sequences associated with all nearest 
30 neighbours of each signal symbol only differ from the bit sequence of that 
signal symbol at one bit position, the enror rate of the transmission system is 
reduced. This form of mapping is referred to as Gray mapping. 

The Invention further relates to an arrangement for generating a reliability 
35 value for a received multilevel signal In relation to a number of predetemiined 
signal symbols each associated with a corresponding bit sequence Including 
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a first bit position; the reliability value being Indicative of likelihood information 
of receiving said multilevel signal; characterised In that the arrangement 
comprises 

- first processing means adapted to identify a first one of the number of 
5 signal symbols as being closest to the received multilevel signal; 

- storage means adapted to store information related to the first signal 
symbol and a second one of the number of signal symbols being closest 
to the first signal symbol and corresponding to a different binary value at 
the first bit position of the respective associated bit sequence than the first 

10 signal symbol; and 

- second processing means adapted to estimate the reliability value on the 
basis of the stored Information. 

The arrangement may be implemented by any processing unit, e.g. a 
15 programmable microprocessor, an application-specific integrated circuit, or 
another integrated circuit, a smart card, or the lil<e. The term processing 
means comprises general- or special-purpose programmable 
microprocessors, Digital Signal Processors (DSP), Application Specific 
Integrated Circuits (ASIC), Programmable Logic Arrays (PLA), Field 
20 Programmable Gate Arrays (FPGA), etc., or a combination thereof. The 
processing means may be a CPU of a computer, a microprocessor, a smart 
card, a SIM card, or the like. The first and second processing means may be 
separate processing means, e.g. separate circuits, or they may be combined 
in one processing means, e.g. performed by suitable Instructions executed by 
25 a programmable microprocessor. 

The term storage means includes magnetic tape, optical disc, digital video 
disk (DVD), compact disc (CD or CD-ROM), mini-disc, hard disk, floppy disk, 
ferro-electric memory, electrically erasable programmable read only memory 
(EEPROM), flash memory, EPROM, read only memory (ROM), static random 
30 access memory (SRAM), dynamic random access memory (DRAM), 
synchronous dynamic random access memory (SDRAM), ferromagnetic 
memory, optical storage, charge coupled devices, smart cards, etc. 



Furthermore, the above discussed features and steps of the method 
according to the invention may be incorporated in the above an-angement 
according to the invention. 

6 The invention further relates to a device for receiving multilevel signals 
comprising an an-angement as described above and in the following. 

The device may be any electronic equipment or part of such electronic 
equipment, >A/here the term electronic equipment includes computers, such 
10 as stationary and portable PCs. stationary and portable radio 
communications equipment. The term portable radio communications 
equipment Includes mobile radio terminals such as mobile telephones, 
pagers, communicators, e.g. electronic organisers, smart phones, PDAs, or 
the like. 

15 

The invention will be explained more fully below in connection with preferred 
embodiments and with reference to the drawings, in which: 

figs, la-b schematically show a receiver according to an embodiment of the 
20 Invention; 

fig. 2 shows an example of a signal constellation with 16 signal symbols; 

fig. 3 shows a flow diagram of a method according to a first embodiment of 
25 the Invention; 

fig. 4 shows an example of a look-up table according to the embodiment of 
fig. 3; 

30 fig. 5 shows a flow diagram of a method according to a second embodiment 
of the Invention; 

fig. 6 shows an example of a look-up table according to the embodiment of 
fig. 5; 
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fig. 7 shows another example of a signal constellation with 16 signal 
symbols; 

fig. 8 shows a flow diagram of a method according to a third embodiment of 
5 the invention; and 

fig. 9 shows an example of a look-up table according to the embodiment of 
fig. 8. 

10 Fig. la schematically shows a receiver according to an embodiment of the 
invention receiving a radio signal from a transmitter via a communications 
channel. Transmitter 101 is adapted to send a signal s via a noisy channel 
102 to receiver 103. The signal s represents one of a set of M signal points 
Si... Sm in a signal space where each signal point is related to a respective bit 

15 sequence of \og2{M) bits. In the presence of noise In the transmission 
channel 102, the receiver 103 receives a signal r" which deviates firom the 
transmitted signal s. In one embodiment the signal is a Code Division 
Multiple Access (CDMA) signal using a spread spectrum technique. The 
receiver 103 comprises a receiver circuit 107 for transforming the received 

20 spread spectrum signal into the signal symbol r. The receiver further 
comprises a channel decoder 108 for decoding the received signal symbol r. 
e.g. a BCJR or Viterbi decoder. The decoder 106 requires soft bit values as 
an input. Hence, the receiver 103 further comprises a circuit 104 which is 
adapted to calculate soft values for the log2(M) bits of the received signal 

25 symbol r and to provide the calculated soft values to the decoder 106. 
According to the invention, the receiver 103 further comprises a memory 105, 
such as on-chip memory, EPROM, flash memory, or the like, in which a look- 
up table is stored for use In an efficient calculation of the soft values by the 
circuit 104, preferably as described in connection with figs. 3-9. 

30 

Fig. 1 b schematically shows a more detailed block diagram of the receiver 
circuit 103 of fig. The circuit 107 comprises a RAKE receiver 110 suitable for 
receiving CDMA signals, i.e. a receiver which uses several baseband 
correlators to individually process several signal multipath components. The 
35 conrelator outputs are combined to achieve improved communications 
reliability and perfonnance (see e.g. "Digital Communications" 4th Edition, by 
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John G. Proakis. McGraw-Hill, 2000). The sampled received radio signal r" is 
fed to the RAKE receiver 110 which generates the signal symbol r to be 
decoded. The circuit 107 further comprises a channel estimator 111 and a 
noise estimator 112. e.g. implementing any suitable channel estimation and 

5 noise estimation technique known in the art. The channel estimator receives 
the received radio signal r*, identifies up to N different radio paths or channel 
taps and estimates corresponding delays Ak. k=1,...,N, and complex channel 
estimate hr = (hri....hrN) of these paths. The channel estimator 111 further 
provides a set of complex combiner weights w = [Wi W2 ... wnI^ to be used by 

10 the rake receiver. Here, denotes a transposed vector. For example, the 
weights may be determined according to an optimisation criterion, such as 
maximising the received signal energy. 

The calculated delays Ak and the combiner weights are provided to the rake 
15 receiver 110. The RAKE receiver 110 comprises delay circuits 115 which 
delay the incoming signal according to the N channel taps. Further, the 
receiver 1 1 0 comprises circuitry 1 1 6 for multiplying the N delayed versions of 
the received signal with a spreading code c for dispreading the spread 
spectaim signals and circuitry 117 for summing the signals to form a radio 
20 symbol. Furthermore, the rake receiver 1 1 0 comprises multiplier circuitry 1 1 8 
for multiplying each of the N radio symbols with the combiner weights (wk)*. 

k=i N, where ()* denotes complex conjugation. Finally, the RAKE receiver 

110 comprises an adding circuit 119 which combines the weighted symbols 
to form the received symbol estimate r which is fed to the soft value 
25 calculation circuit 1 04. 

When using a multilevel constellation of signal points Si...Sm the amplitude 
information should be maintained in order to ensure successful demodulation 
in the receiver. Consequently, the reference points Si...Sm should be scaled 

30 properly. In the following it is assumed that the channel estimator 111 
estimates the channel gain on the basis of a reference channel hr which has 
a channel gain that may be different from the actual gain of the traffic 
channel, e.g. a HS-DSCH. The gain difference between the reference 
channel and the traffic channel may be denoted with g. Hence, the received 

35 symbol r after the RAKE receiver 1 1 0 may be expressed as 
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r=g w"hrS + n, 

where w" is the Hermltian conjugate of w, w"hr denotes an Inner product, s 
is the transmitted symbol and n Is a noise term, e.g. representing additive 
5 white Gaussian noise (AWGN). The gain parameter g is signalled to the 
receiver, w is selected by the combiner in the receiver, and hr are the 
channel estimates. Hence, at the receiver, the reference signal symbols 
Si...Sm may be scaled appropriately, according to 

10 $j=gw"hrSj, j=1.....M. (0) 

In fig. 1b. the receiver circuit 107 comprises circuit 1 13 adapted to calculate 
the above scaling factor g w"hr and a multiplier circuit 1 14 for multiplying the 
reference symbols with the scaling factor, resulting In properly scaled signal 
15 symbols Si,...,S|y/|. which are fed into the soft value calculation circuit 104. 

Finally, the noise estimator 1 12 provides an estimate of the signal noise level 
a which Is fed Into the soft value calculation circuit 104. 

20 According to the Invention, the soft value calculation circuit locates the signal 
symbol which is closest to the received signal r and calculates corresponding 

soft values Lm for bit m, m=1 log2(M), e.g. according to one of the 

embodiments discussed in connection with figures 3-9. 

25 It Is noted that the receiver circuit described in connection with figs, la-b 
merely serves as an example, and the scope of the Invention Is not limited to 
the type of receiver, nor to the above scaling of signal symbols. 

Fig. 2 shows a signal constellation with 16 signal symbols. The signal 
30 constellation comprises M=16 signal points Si through Sie in a two- 
dimensional signal space, e.g. the l/Q components in a 16QAM signal 
constellation. Preferably, the signal points are distributed regularly, such that 
the distance to the nearest neighbours of each signal point is the same. The 
reference points may take values that suit the Implementation in question. 
35 However, alternatively, other signal constellations may be chosen. In fig. 2, 
16 different bit sequences 0000 through 1111, each consisting of log2(M)=4 
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blts, are mapped onto the signal points Si-Sie. Preferably, tine mapping of tine 
bit sequences to the signal points is chosen such that the bit sequence of 
each signal point only differs firom those of the nearest neighbours by one bit, 
thereby optimising the decoding performance. For example, in fig. 2, signal 
5 point Ss has three nearest neighbours. S4. S7, and S12. The bit sequence of 
S4, i.e. 1 1 1 1 . differs from the sequence 1 1 1 0 of Ss only at bit position 4, etc. 
Altennatively. other mappings may be chosen. 

For every bit position m mapped on a signal point, the signal points in the 
1 0 constellation may be divided into two sets where the signal points in each set 
have bit value 0 and 1, respectively, at that position. In the following, the set 
of signal points with a 0 at the m-th position Is denoted Ao.m. and the 
con-esponding set with a 1 at the m-th position is denoted Ai.ni. For example, 
In the example of fig. 2 and for m=1, Ai.iKSa, S4, S7, Ss. Sn, S12, S15. Sie} 
16 and Ao.i={Si. S2. Ss. Se. S9. S10. S13. S14}. The sets are of equal size with 
M/2 elements each. 

When one of the signals Si...., Sie is transmitted over a noisy channel, the 
received signal will differ from the transmitted signal according to a 
20 corresponding distribution. The actual shape and width of the distribution of 
received signals depends on the characteristics of the noise. In fig. 2, the 
cross 201 represents an example of a received signal r. 

Prior to providing the received 16QAM radio symbols to a decoder, e.g. a 
25 turbo decoder, they are converted into soft values. Hence, a soft value is 
calculated for each bit of every 16QAM symbol. A soft value of the m-th bit in 
the sequence mapped to r may be defined as 

. P(sm=t) . P(sm=l)P(r|^m=l) _.^,gHfm^ (D 

where Sm Is the m-th bit in the bit sequence represented by the transmitted 
signal, and P(Sm=i| r), 1=0.1, are the a posteriori probabilities of the bit Sm 
where r Is the received signal. It is noted that the second equality assumes 
35 that Sm=1 and s,„=0 are equally probable in the chosen alphabet. Otherwise, 
the overall ratio of probabilities should be taken into consideration in the 
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10 



15 



following. However, this would only give rise to a constant factor. Hence, Lm 
corresponds to a log-|j|<elihood ratio of probabilities. The probabilities 
P(r|Sm=i) in eqn. (1 ) may be written as 

P('-|Sm=0=^ ZKr.s) 1 = 0.1. (2) 

Hence, tlie calculation of tlie above probability involves a summation over 
M/2 terms eachi including a joint probability P(r,s), This is a computationally 
expensive task, especially if M is large, e.g. M=64. 

In many applications, the above soft values Ui may be approximated by 

max P(r.s) ,3. 
max P(r.s) ''p(r|§o.m) 

seA 



0,m 



where §i,m, i=0,1, are the signal points that result in the largest contribution to 
the sums In eqn. (2). Hence, In the calculation of the probabilities, the sums 
over M/2 terms are approximated by the their respective dominant terms. 
20 according to 



'og 2P(r.s)«log max P(r,s) 

The above approximation is often referred to as the "max log MAP" 
25 approximation which yields a good approximation in cases where the above 
sums are dominated by one term, as for example in the case of Gaussian 
noise when the signal to noise ratio (SNR) is large. The above probabilities 
depend on the distances between the received signal r and the respective 
signal points. For example, in the case of additive zero-mean Gaussian noise 
30 with variance o^, the log-likelihood ratio of eqn. (3) may be expanded as 



0-2 



Lm=iog— ■ 

a ^exp 



exp(^-|r-Si,m|2/a2j 



36 



r7~T~1^7^ = ^ U'-^0''"l -|''-Hm| ) (5) 
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It is noted that it Is assumed that the §i.ni are scaled corresponding to the 
received signal according to equation (0) above. In the following, we define 
di.m=|r- §i.m|. for i=0.1, to be the distances between the received signal r and 
the closest signal points in the sets Ai.m, respectively. For example, in fig. 2, 

5 for m=1, do,i corresponds to the distance 82 between r and the closest signal 
point in A0.1, i.e. Se, while di.i corresponds to the distance 81 between r and 
the closest signal point in A1.1, i.e. Ss. According to the invention, the 
likelihood ratio in equation (5) is obtained by first Identifying the closest 
signal point Se, and then determining the distances 5i and 82, as will be 

10 described in greater detail in connection with figs. 3-4. Hence, according to 
the invention, a computationally expensive calculation of all the distances 
between r and all the signal points S1...S16 in order to identify the shortest 
distances 81 and 82 is avoided. Alternatively to identifying Se by means of a 
look-up table and then calculating 82, the distance 83 between Se and Se may 

15 be looked up and used as an approximation instead of 82, thereby saving 
additional computational resources. This will be described in greater detail in 
connection with figs. 5-6. A further embodiment of the invention will be 
described In connection with figs. 8-9. 

20 It is noted that, preferably, in the above estimation of the reliability values, a 
proper scaling of the signal points in the QAM constellation is taken into 
consideration. If this scaling is taken into consideration, the above log- 
likelihood ratio may be written as 

2S L„=a-2.(lr-So.mr-hHmr} 

As described in connection with fig. 1b, when using a multilevel constellation 
as in the example of fig. 2, the amplrtude information should be maintained In 
30 order to ensure successful demodulation in the receiver. Consequently, the 

reference points Sj, j=1 16 should be scaled properly. If this scaling Is 

taken into consideration, the above log-likelihood ratio may be written as 

Lm=K{|r-So,mf-|r-§,.„,f) (7) 

35 

i.e. with the properly scaled signals 
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Si.m=— TT- m = 1...Jog2(M).i = 0,t (8) 
gw"hr 



5 and where 

K = ^w"hrf I<P- (9) 
is a constant which depends on the signal to noise ratio. 

10 

Fig. 3 shows a flow diagram of a method according to a first embodiment of 
the invention. According to this embodiment, the soft values Lm are 
calculated using the approximation in equation (7). Initially, in step 301, a 
signal r is received and, in step 302, the signal point s from the set of signal 

15 points Si...Sm which, according to a Euclidean metric, is closest to r is 
identified. For example, an efficient way of identifying § is by means of a 
slicer. In step 303, the distance 5i between r and § is calculated. 

Subsequently, for bit positions m-1 log2(M), the following steps are 

performed: In step 304, the signal point s which is closest to s is looked up in 

20 a look-up table 308. This signal point corresponds to the signal point which is 
closest to r and has the opposite bit value at position m than §. In step 305, 
the distance 52 between r and § is calculated. Based on the distances 5i and 
52, the soft value Lr, is now approximated according to eqn. (7) above: If the 
bit value §m of § at position m is 0, the soft value is approximated by Lm- 

25 K (52)^-(5i)2 (step 306). Othenvise the soft value is approximated by Lm= 
K (5i)2-(52)^ (step 307). Here, K Is a constant which depends on the noise 
distribution as described above. Referring to the example illustrated in fig. 2, 
the closest signal point to the received signal r (marked by the cross 201) is 
s^Ss. When calculating a soft value Li for the first bit position m-1 using the 

30 method of fig. 3, the first bit in Ss is identified to be §i = 1. From a pre- 
computed look-up table, e.g. as illustrated in fig. 4, the closest signal point 
with a "0" In the first bit position is S-Se. Hence, the distances 5i and 52 may 
be calculated as 5i=:|r-S8i and S2=|r-S6|, respectively, where |-| denotes the 
Euclidean distance. Thus, the soft value Li is approximated by Li = K'(52)2- 

35 (5i)2 
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Consequently, the method according to this embodiment requires at the most 
1 + Iog2(M) distance calculations, since the closest distance 5i has to be 
calculated once for a received signal (step 303) and. for each bit position m, 
the distance 62 is calculated in step 305. This is to be compared with M 
5 distance calculations when all distances to all signal points are calculated. 
Hence, the computational complexity of this method only grows 
logarithmically with the size of the symbol alphabet rather than proportional to 
the alphabet size. This is a considerable reduction of the computational 
complexify, in particular for large alphabet sizes. 



Fig. 4 shows an example of a look-up table according to the embodiment of 
fig. 3. The look-up table 308 identifies the pre-computed closest signal points 
and is indexed by the bit numbers m and the signal points Si...Sm. Each row 
corresponds to one of the signal points Si...Sm. For example, row 402 

1 5 corresponds to signal point S2, such that each element In row 402 identifies a 
signal point which is the closest to S2 among all signal points having a bit 
value opposite to S2 at bit position m. Each entry in table 308 consumes 
log2(M) bits for identifying one out of M signal points. Furthemriore, the table 
consists of M rows and log2(M) columns. Consequently, the table requires M' 

20 • [\oQ3iM)f bits. For example for M=8 the memory consumption Is 72 bits and 
for M=16 the memory consumption is 256 bits. Hence, it is a further 
advantage of this embodiment that it only requires little storage space for 
storing the look-up table. It is a further advantage of this embodiment, that It 
yields a good approximation of the soft values, thereby providing a good 

25 decoding performance. 

Fig. 5 shows a flow diagram of a method according to a second embodiment 
of the Invention. Again, this embodiment utilises the approximation of 
equation (7) for the calculation of the soft values Lm. As in the embodiment of 

30 fig. 3, in the initial step 501 . a signal r is received and, In step 502. the signal 
point S from the set of signal points Si...Sm, which is closest to r is identified, 
e.g. by means of a slicer. In step 503, the distance Si between r and s Is 
calculated. Subsequently, for bit positions m=1,...,log2(M), the following 
steps are performed: In step 504, the distance S3 between § and the signal 

35 point s which is closest to s and has the opposite bit value at position m is 
looked up in a look-up table 508. Subsequently, this distance 83 Is used as an 



10 
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approximation for tine distance 82 between r and s when approximating the 
soft value Lm according to eqn. (7) above. Hence, if the bit value §m of § at 
position m is 0, the soft value Is approximated by Lni= K (53)^-(5i)2 (step 
506). Otherwise the soft value is approximated by Lm= K-(6i)2-(53)2 (step 
5 507). Again, K is a constant which depends on the noise distribution. 
Referring again to the example Illustrated in fig. 2, the closest signal point to 
the received signal r is §=S8. When calculating a soft value Li for the first bit 
position m=1 using the method of fig, 5, the first bit in Ss is Identified to be §1 
= 1. From a pre-computed look-up table, e.g. as illustrated In fig. 6, the 
1 0 distance to the closest signal point with a "0" in the first bit position Is di.8=53. 
Hence, the distance 5i is calculated as 5i=|r-S8| and 82 is approximated by 63. 
Thus, the soft value Li is approximated by Li = K (83)2-(8i)2. 

Consequently, as the pre-computed distance 83 is used as an approximation 
16 for 82, the method according to this embodiment requires only one distance 
calculation, i.e. the calculation of 81 (step 503). Again, this Is to be compared 
with M distance calculations when all distances to all signal points are 
calculated. 

20 Hence, it Is an advantage of this embodiment that the computational 
complexity does not grow with the size of the symbol alphabet, thereby 
yielding a computationally efficient method of approdmating soft values. 

Fig. 6 shows an example of a look-up table according to the embodiment of 
26 fig. 6. The look-up table 608 comprises the pre-computed distances between 
the signal points in a constellation of size M. A distance dm,k in table 508 
denotes the Euclidean distance between signal point Sk and the closest 
signal point with opposite bit value at position m. Assuming that each 
distance is stored with a resolution requiring ti bits, each entry in the table 
30 608 requires -n bits. Furthermore, the table consists of M rows and log2(M) 
columns. Consequently, the table requires it M • log2(M) bits. For example, 
for M=8 the memory consumption is 241^ bits and for M=16 the memory 
consumption Is 64ti bits. Hence, it Is a further advantage of this embodiment 
that It requires little storage capacity. In an emtKJdIment where the resolution 
36 of the pre-computed distances is higher than log2(M) bits. I.e. T|>log2(M), 
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processing time is traded for memory space in comparison with the 
embodiment of figs. 3-4. 

it is noted that additional storage space may be saved by only storing each 
5 distance once, i.e. in case the same distance appears in two or more entries 
of the table, a reference to that distance may be stored In one of the entries, 
Instead. 

Alternatively, other layouts of a look-up table may be used. For example, in 
10 one embodiment, the look-up table may comprise all M(M-1)/2 mutual 
distances between the signal points Sk. thus requiring 'nM(M-1)/2 bits of 
storage. However, for M>4 this embodiment requires larger storage capacity 
than the embodiment of fig. 6. 

16 Fig. 7 shows another example of a signal constellation with 16 signal 
symbols. As In fig. 2, the signal constellation comprises M=16 signal points 
Si through Sie in a two-dimensional signal space, e.g. the l/Q components In 
a 16QAM signal constellation. The signal points are distributed regularly, 
such that the distance to the nearest neighbours of each signal point is the 

20 same. In the example of fig. 7. they are assumed to be selected such that 

Sk = Xk + j yk, where Xk. yk e [-3d, -d, d, 3dJ, k-1,...,M, 

where d is an arbitrary constant and where f = -1 . For example, d may be 
25 chosen to d=:1. However, alternatively, other signal constellations may be 
diosen. 

In fig. 7, 16 different bit sequences 0000 through 1111, each consisting of 

log2(16)=4 bits, are mapped onto the signal points Si Sie. Preferably, the 

30 mapping of the bit sequences to the signal points Is chosen to be a Gray 
mapping. I.e. such that the bit sequence of each signal point only differs from 
those of the nearest neighbours by one bit, thereby optimising the decoding 
performance. 
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As above, for every bit position m mapped on a signal point, the signal points 
in the constellation may be divided into two sets Ao.m and Ai.m, where the 
signal points in each set have bit value 0 and 1, respectively, at that position. 

5 Fig. 8 shows a flow diagram of a method according to a third embodiment of 
the Invention. Again, this embodiment utilises the approximation of equation 
(7) for the calculation of the soft values Lm- As In the embodiment of fig. 3, in 
the initial step 801. a signal r is received. The received signal may be written 
as r = Re(r) + j lm(r) and, in the following the magnitude of the I- and Q 

10 components of r will be denoted by a = | Re(r) | and b = | lm(r) |. respectively. 
After the received symbol is combined in the combiner, in step 802, the 
signal point s from the set of signal points Si...Sm, which is closest to r is 
Identified. In this embodiment it is assumed that the constellation of signal 
points corresponds to the constellation of fig. 7. In fig. 7, each signal point 

15 corresponds to a decision region where the decision regions are separated 
by a set of decision boundaries 701 through 706. Hence, the closest signal 
point s may be found by performing two comparisons of the inphase 
component and the quadrature component, respectively. For example, if 
Re(r) < 0 (decision boundary 705) and Re(r) <-2d (decision boundary 706) 

20 and If lm(r)>0 (decision boundary 702) and lm(r)>2d (decision boundary 
701), the received signal lies in the decision region corresponding to Si, i.e. 
§ = Si Is the closest signal point. Subsequently, for each bit positions 
m=1,...,log2(16)=1,...,4, the soft value U may be calculated using the 
approximation of eqn. (7), assuming proper scaling. Consequently, In this 

25 example the soft value Li for the first bit Is 

Li(Si) = K(|-a+j b - (d+3j d)p - |-a+j b - (-3d+3j d)|2) ^ ^ (8ad-8d2). 

Hence, in the above equation, instead of computing two distances squared, 
30 each Involving a calculation of the type |x+jyp, the soft value may be 
calculated by scaling the inphase amplitude a of the received symbol with 
8dK and, subsequently, by adding a constant -8Kd^ It Is further noted that 
the constant d may be chosen as any suitable positive real number. 

35 The remaining three soft values for a received signal In the decision region 
corresponding to Si are accordingly 




L2(Si) = K(|-a+jb - (-3cl+3jd)p - |-a+jb - (-3cl-jd)|^) = K (-Sbd+Sd^) 
USi) = K(|-a+j b - (-d+3j d)|2 _ ^^^^y^ - (.3d+3j d)p) = K (4ad-8d2) 
L4(Si) = K(|-a+Jb - (-3d+3jd)|2 - |-a+jb - (-3d+jd)|2) = K (4bd-8d2), 

5 as in the constellation of fig. 7 the closest symbols with opposite second, 
third, and fourth bit compared to Si are Sg = -3d - jd, S2 = -d + 3jd, and 
Ss = -3d + jd, respectively. 

The table 808 of fig. 9 illustrates the calculated soft values for all decision 

10 regions corresponding to the symbols Si,...,Si6, and for all bits, m=1,...,4. As 
can be seen from table 808, all soft values may be calculated by scaling one 
of the inphase component a or quadtrature component b of the received 
signal r and subsequently adding a constant. Hence, using the pre-calculated 
equations of table 808, the soft values may be calculated in a very efficient 

16 way. In one implementation, each entry of the look-up table 808 may 
comprise the scaling factor, the constant to be added and a bit indicating 
whether it is the inphase component a or the quadtrature component b of the 
received signal r which is to be scaled for a given soft value. Preferably, the 
table is indexed by the decision region and the bit values. It is noted, 

20 however, that many of the entries of table 808 are identical. Consequently, it 
will be apparent to a skilled person that table 808 may be stored in a memory 
efficient manner, e.g. by storing a list of the distinct entries and, in table 808. 
referring to the corresponding list members. In general, it is noted that 
constellations which are Gray coded or show another regularity, the 

25 redundancy of the entries in table 808 may be utilised to reduce the memory 
consumption of table 808. 

Refening again to fig. 8, in steps 804-805 the soft values for the identified 
decision region and for all bits are calculated. In step 804, the relation to be 
30 calculated, i.e. the scaling factor and the constant to be added, are retrieved 
firom a stored table 808 in memory, e.g. a look-up table as shown in fig. 9. 
The retrieved relation is calculated in step 805 resulting in the soft value for 
the corresponding bit number. 

35 It is noted that the processing load in the receiver may further be decreased 
by pre-calculating the relations of table 808 and by storing the pre-calculated 
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soft values: Assuming that the inphase and the quadrature components of 
the received signal each are quantised to n bits, the soft values of table 808 
may be precalculated and tabulated for every different inphase and 
quadrature value, thereby further decreasing the required calculations, as the 

5 scaling and adding of step 805 are not necessary in this embodiment. 
However, such a table of pre-calculated soft values Increases the memory 
consumption. Above, a=|Re(r)| and b=|lm(r)| were defined as the absolute 
values of the real and imaginary parts of r, respectively, i.e. without sign 
information. Hence, a and b, each are represented by n-1 bits. If each of the 

10 pre-calculated soft values Is to be represented by m bits, the total memory 
consumption of a full table is m-2"'^-4-16 bits (for each of the 16 decision 
regions and each of the 4 bits, 2""^ different soft values are stored, each with 
a precision of m bits). For example, for n=4, the total memory is 512m bits. 
Note, however, that in this embodiment, the pre-calculated table still needs to 

1 5 be multiplied with the factor K 

It Is noted, that the above memory consumption may be further reduced by 
utilising the fact that many of the entries of table 808 are Identical and by 
utilising the symmetry of the constellation of fig, 7. 

20 

It Is noted that the invention was described in connection with soft values 
defined as a log-llkelihood ratio Indicating a reliability value for the bit values 
of a received sequence. However, other definitions of soft values depending 
on the distance of the received signal to the signal points may be used as 
25 well. 

It Is further noted that the signal constellations of figs. 2 and 7 are merely 
used as examples. The calculation of soft values according to the Invention Is 
not limited to these signal constellations. 

30 

Finally, It is noted that the embodiments described in connection with figs. 3- 
6 are particularly well suited for large signal constellations, as they save 
memory, whereas the embodiment of figs. 8-9 is particularly well suited for 
implementing a medium-size signal constellation, e.g., 16QAM, as the save 
35 computational resources. 
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CLAIMS 



1. A method of generating a reliability value {Lk,m) for a received multilevel 
signal (r) in relation to a number of predetermined signal symbols (Si,....Sm) 

5 each associated with a conresponding bit sequence including a first bit 
position (m); the reliability value being indicative of likelihood infomnation of 
receiving said multilevel signal; 

characterised in that the method comprises the steps of 

- identifying (302; 602; 802) a first one (§k) of the number of signal symbols 
10 as being closest to the received multilevel signal; and 

- estimating (306. 306, 307; 606, 607; 804, 806) the reliability value based 
on stored information related to a second one (sr) of the number of signal 
symbols being closest to the first signal symbol and corresponding to a 
different binary value at the first bit position of the respective associated 

1 5 bit sequence than the first signal symbol. 

2. A method according to claim 1, characterised In that the stored infomnation . 
comprises an identification of the second signal symbol, and the step of 
estimating the reliability value further comprises the steps of 

20 - determining (303) a first distance (5i) between the received signal and the 

first signal symbol; and 

- determining (306) a second distance (62) between the received signal and 
the second signal symbol. 

26 3. A method according to claim 1 . characterised In that the stored infomnation 
comprises a second distance (83) between the first signal symbol and the 
second signal symbol, and the step of estimating the reliability value further 
comprises the step of determining (603) a first distance (81) between the 
received signal and the first signal symbol. 

30 

4. A method according to claim 2 or 3, characterised In that the step of 
estimating the reliability value comprises the step of determining (306. 307; 
606, 607) a polynomial function of the first distance and the second distance, 
multiplied by a predetenmined constant (K). 

36 
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6. A method according to claim 4, characterised in that the predetermined 
constant Is selected depending on the noise distribution of the received 
multilevel signal. 

5 6. A method according to claim 1 , characterised in that the stored information 
is indicative of a number of functional relations between the received 
multilevel signal and the reliability value, and the step of estimating the 
reliability value further comprises the step of selecting 804 a functional 
relation of said number of functional relations dependant on the first signal 
1 0 symbol and the first bit position. 

7. A method according to claim 1. characterised in that the stored 
information comprises for each signal symbol and bit position an 
approximation of the corresponding reliability value. 

16 

8. A method according to any one of the claims 1 through 7, characterised in 
tiiat the stored information is stored in a look-up table (308; 508; 808) 
indexed by the number of signal symbols and the bit positions of the number 
of bit sequences. 

20 

9. A method according to any one of the claims 1 through 8, characterised in 
that the method further comprises the step of providing the reliability value as 
an input to a decoder (1 06). 

25 10. A method according to any one of the claims 1 through 9, characterised 
in that the first signal symbol is identified by means of a slicer. 

11. A method according to any one ofVne claims 1 through 10, characterised 
in that the likelihood information comprises a log-likelihood ratio. 

30 

12. A method according to any one of tiie claims 1 through 11, characterised 
in that the step of identifying the first signal symbol as being closest to the 
received multilevel signal comprises tine step of Identifying the first signal 
symbol as being closest to the received multilevel signal with respect to a 

35 Euclidean distance measure in a signal space. 
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13. A method according to claim 12, characterised in that the signal space is 
related to the complex plane in quadrature amplitude modulation. 

14. A method according to any one of the claims 1 through 13, characterised 
5 in that the number of signal symbols are associated with the number of bit 

sequences such that the bit sequences associated with all nearest 
neighbours of each signal symbol only differ from the bit sequence of that 
signal symbol at one bit position. 

10 15. An arrangement for generating a reliability value (Um) for a received 
multilevel signal (r) in relation to a number of predetermined signal symbols 
(Si,...,Sm) each associated with a corresponding bit sequence including a 
first bit position (m); the reliability value being indicative of likelihood 
Information of receiving said multilevel signal; 

15 characterised in that the arrangement comprises 

- first processing means (1 04) adapted to Identify a first one of the number 
of signal symbols as being closest to the received multilevel signal; 

- storage means (105) adapted to store information related to the first 
signal symbol and a second one of the number of signal symbols being 

20 closest to the first signal symbol and corresponding to a different binary 
value at the first bit position of the respective associated bit sequence 
than the first signal symbol; and 

- second processing means (104) adapted to estimate the reliability value 
on the basis of the stored information. 

25 

16. An arrangement according to claim 15, characterised in that the storage 
means Is adapted to store an Identification of the second signal symbol, and 
the second processing means is ftjrther adapted to determine a first distance 
between the received signal and the first signal symbol and a second 

30 distance between the received signal and the second signal symbol. 

17. An arrangement according to claim 15, characterised in that the storage 
means is adapted to store a second distance between the first signal symbol 
and the second signal symbol, and the second processing means is further 

35 adapted to determine a first distance between the received signal and the 
first signal symbol. 
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18. An arrangement according to claim 16 or 17, characterised in that tlie 
second processing means is furtlier adapted to detemnine a polynomial 
function of the first distance and the second distance, multiplied by a 
predetermined constant. 

6 

19. An anrangement according to claim 18, characterised in that the 
predetermined constant is selected depending on the noise distribution of the 
received multilevel signal. 

10 20. An anrangement according to claim 15, characterised in that the storage 
means is adapted to store information indicative of a number of functional 
relations between the received multilevel signal and the reliability value, and 
the second processing means is further adapted to select a functional 
relation of said number of functional relations dependant on the first signal 

1 5 symbol and the first bit position. 

21. An arrangement according to claim 15, characterised in that the stored 
information comprises for each signal symbol and bit position an 
approximation of the corresponding reliability value. 

20 

22. An arrangement according to'any one of the claims 15 through 21, 
characterised in that the storage means is adapted to store the stored 
information as a look-up table indexed by the number of signal symbols and 
the bit positions of the number of bit sequences. 

25 

23. An arrangement according to any one of the claims 15 through 22, 
characterised in that the first processing means further comprises a slicer. 

24. An arrangement according to any one of the claims 15 through 23, 
30 characterised in that the likelihood information comprises a log-likelihood 

ratio. 

25. An arrangement according to any one of the claims 15 through 24, 
characterised in that the first processing means is further adapted to identify 

35 the first signal symbol as being closest to the received multilevel signal with 
respect to a Euclidean distances in a signal space. 
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26. An arrangement according to claim 26, characterised In that the signal 
space Is related to the complex plane in quadrature amplitude modulation. 

5 27. An arrangement according to any one of the claims 16 through 26, 
characterised in that the number of signal symbols are associated with the 
number of bit sequences such that the bit sequences associated with all 
nearest neighbours of each signal symbol only differ from the bit sequence of 
that signal symbol at one bit position. 

10 

28. A device (103) for receiving multilevel signals comprising an arrangement 
according to any one of the claims 1 5 through 27, 

29. A device according to claim 28, characterised in that the device further 
15 comprises a decoder (106) adapted to receive an input signal from the 

an-angement indicative of the detemnined reliability value. 



30. A device according to claim 28 or 29, characterised in that the device Is a 
mobile temnlnal. 
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Soft value calculation for multilevel signals 



Abstract 
5 

Disclosed is a method of generating a reliability value for a received 
multilevel signal (r) in relation to a number of predetermined signal symbols 
(Si,...,Si6) each associated with a corresponding bit sequence Including a 

10 first bit position; the reliability value being indicative of likelihood information 
of receiving said multilevel signal. The method comprises the steps of 
identifying a first one (Ss) of the number of signal symbols as being closest to 
the multilevel signal; and estimating the reliability value based on stored 
information related to a second one (Se) of the number of signal symbols 

15 being closest to the first signal symbol and corresponding to a different binary 
value at the first bit position of the respective associated bit sequence than 
the first signal symbol.. 

20 (Figure 2 should be published). ... 
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Defective images within this document are accurate representations of the original 
documents submitted by the applicant. 



Defects in the images include but are not limited to the items checked: 



J-M^ACK BORDERS 

□ IMAGE CUT OFF AT TOP, BOTTOM OR SIDES 

□ FADED TEXT OR DRAWING 

□ BLURRED OR ILLEGIBLE TEXT OR DRAWING 

□ SKEWED/SLANTED IMAGES 

□ COLOR OR BLACK AND WHITE PHOTOGRAPHS 

□ GRAY SCALE DOCUMENTS 

□ LINES OR MARKS ON ORIGINAL DOCUMENT 

□ REFERENCE(S) OR EXHIBIT(S) SUBMITTED ARE POOR QUALITY 

□ OTHER: 

IMAGES ARE BEST AVAILABLE COPY. 
As rescanning these documents will not correct the image 
problems checked, please do not report these problems to 
the IFW Image Problem Mailbox. 
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